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A gamma-ray burst fireball is likely to contain an admix- 
ture of neutrons, in addition to protons, in essentially all pro- 
genitor scenarios. Inelastic collisions between differentially 
streaming protons and neutrons in the fireball produce v^(u^) 
of ^10 GeV as well as i^e(j^e) of ~ 5 GeV, which could produce 
~ 7 events/year in km^ detectors, if the neutron abundance 
is comparable to that of protons. Photons of ~ 10 GeV from 
tt" decay and ~ 100 MeV Ve from neutron decay are also 
produced, but will be difhcult to detect. Photons with ener- 
gies ^ 1 MeV from shocks following neutron decay produce 
a characteristic signal which may be distinguishable from the 
proton-related MeV photons. 



I. INTRODUCTION 

Gamma-ray burst (GRB) sources are distributed 
throughout the universe and their energy output is mea- 
sured to he. a substantial fraction of a solar rest mass 
equivalent A variety of observations support the in- 
terpretation that these events are caused by cataclysmic 
stellar collapse or compact mergers, producing a fireball 
with bulk expansion Lorentz factor P ~ 10^ - 10^. In 
the standard GRB model a fireball made up of 7, e^ and 
magnetic fields with an admixture of baryons is produced 
by the release of a large amount oLenergy E <; 10^'^ ergs 
in a region Tq ^ \Q^ro7 cm (e.g.EI). The observations 
indicate that typical fireballs are characterized by a lu- 
minosity L ~ 10^^L52 erg s~^ and durations tw = lOt^ji s 
in the observer frame, with a large spread in both quanti- 
ties. The outflow is controlled by the value of the dimen- 
sionless entropy rj — (LjMc?) injected at Tq- Previous 
discussions of fireball models have generally focused on 
the charged particle components, since they determine 
directly the photon signal. However, consideration of a 
neutron component introduces qualitatively new effects. El 
In a p, n fireball, for values of 77 ^ 400, the neutrons and 
protons acquire a relative drift velocity causing inelastic 
71, p collisions and creating neutrinos. 

We investigate here the neutrino and photon signals 
from n,p collisions following decoupling in GRB. The 
~ 10 GeV neutrinos from this mechanism depend upon 
the presence of neutrons in the original explosion, but the 
neutrinos are created in simple physical processes occur- 
ring in the later stages of the fireball. Oii-the other hand, 
the 10'"' GeV neutrinos discussed in refs.a require the ac- 



celeration in shock waves of ultra-high energy protons 
interacting with photons. Thus the 10 GeV and the 10^ 
GeV neutrinos reflect very different astrophysical pro- 
cesses and uncertainties. Other processes, e.g. neutrinos 
from p,p collisiong3 also require shocks but have lower 
efficiencies, while 10-30 MeV neutrinos from the original 
explosiodj are much harder to detect due to the lower 
cross sections. 



We show (§111) that the 10 GeV neutrinos could be 
detectable by future km^ size detectors. The associated 
^ 10 GeV 7-ray fluences are compatible with current de- 
tection rates, and may be detectable with future space 
missions. The dependence of these signals on the neu- 
tron/proton ratio ^ provides a new tool to investigate 
the nature of the GRB progenitor systems. Moreover 
the predicted neutrino event rate depends on the asymp- 
totic bulk Lorentz factor of the neutrons, which is linked 
to that of the protons. The latter affects all the electro- 
magnetic observables from the GRB fireball, including 
the photospheric and shock radii, as well as the particle 
acceleration and non-thermal photon production. 



II. DYNAMICS, N-P DECOUPLING AND PIONS 

Above the fireball injection radius Tq the outflow ve- 
locity increases through conversion of internal energy 
into kinetic energy, the bulk Lorentz factor P varying 
as P ^ Tq/T" = '"/''o, where T' is the comoving tem- 
perature and T'^ = \. 21^12 oi^"^ MeV is the initial tem- 
perature at To (henceforth denoting with primes quan- 
tities measured in the comoving frame). The flow may 
be considered spherical, which is a valid approximation 
also for a coUimated outflow of opening angle Bj > P~^, 
for the conditions discussed here. In a pure proton out- 
flow the linear growth of P saturates when it reaches an 
asymptotic value P/ < ?7 ~ constant, the value r\ be- 
ing achieved when the fireball converts all its luminosity 
into expansion kinetic energy. For an n,p fireball, be- 
yond the injection radius Tq the comoving temperature 
is low and nuclear reactions are rare, so the njj) ratio ^ 
remains constant. Since the thermal velocities are non- 
relativistic, decoupling of the n and p fiuids is essential 
for high-energy neutrino production. 

At the base of the outflow the n and p components 
are coupled by nuclear elastic scattering. In terms of 
the CM relative energy trei and the relative velocity 

UoC. The CM en- 
"'^^ is approximately 



between nucleons, a'^iv'^ 
ergy dependence oc e' 



el ^rel 
-1/2 , 



el 
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valid between energies MeV and the pion produc- 
tion threshold ~ 140 MeV, and CTq ti^ 3 x 10~26 
cm^ is the pion formation cross section above thresh- 
old. The p and n are cold in the comoving frame, and 
remain well coupled until the comoving n,p scattering 
time t'j^p ~ {n'p<7oc)~^ becomes longer than the comov- 
ing expansion time t'^^^ ^ r/cT. Denoting the comoving 
neutron density njj — ^n'p with ^ ^ 1, mass conservation 
implies Tip = L/[{l+£_)4:nr'^mpC^rri]. The n,p decouphng 
occurs in the coasting or accelerating regimes depending 
on whether the dimensionless entropy rj is below or above 
the critical value 

?77r = {L(7T,/4:T:mpC^ro{l + ^y^^ 

^ 3.9xio^Lli\;,'^U[i + &2)-'/' ■ (1) 

Figure [l] shows the dependence of F on radius for different 
r]. For low values, rj <, rj^^, the condition t'^-^^ ^ t'^^^ is 
achieved at a radius rnp/vo = 'H-niv-n / "hY , which is beyond 
the saturation radius rg/ro ^ rj at which both n and p 
start to coast with F ~ 77 = constant. In this case, even 
after decoupling both n and p continue to coast together 
due to inertia, and their relative velocities never reach 
the threshold for inelastic collisions. 
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FIG. 1. Schematic behavior of the bulk Lorentz factor F 
as a function of radius r for various values of the dimensionless 
entropy rj, the decoupling radius r„p being indicated with a 
diagonal slash. Curve 1 is for rj < rjTr, where the n and p 



achieve the same asymptotic r„/ 



-pf 



rj. Curve 2 is for 



T] > titt, and in this case n,p decoupling occurs before protons 
have reached their asymptotic Lorentz factor, which is larger 
than that of neutrons. This leads to inelastic n,p collisions, 
pion formation and neutrino emission at r„p,2. 

For 77 ^ 777r, on the other hand, the n,p decoupling con- 
dition t'^p ^ t'^^p occurs while the protons (and neutrons) 
are still accelerating as Fp ~ (r/ro), at a radius 



(rnp/ro) = rj-ni'n/'n-n) , for ?7 > 77^ 



Beyond this decoupling radius the p can still continue to 
accelerate with Fp oc r (as long as they remain coupled 



to the photons). However the neutrons are no longer 
accelerated, since they only interact with the protons, 
and they continue to coast with the value of F ~ F„^ ~ 
constant achieved up to that point, 

F„/= (3/4)?7,(77/r?.)-i/3 (fo^ ^ t^^) 

3 X WL]!,\-J'' ([1 + ivM-'/' , (3) 

where the (3/4) factor comes from a numerical solution!! 
of the coupling equations. 

When the n,p decoupling condition r/ ^ ry^r is satisfied, 
the relative n,p drift velocity Vrei c and the inelastic 
pion production threshold e' > 140 MeV is reached. Since 
ctq ~ (T^ , the condition t'^p ^ t'^^^ implies that the optical 
depth to pion formation is of order unity. Thus, for rj ^ 



tItt the radius 

' np 



is not only a decoupling radius 



but also an effective "pionospheric" radius. 

The lowest energy threshold processes at are 



p + n — >p + p + TT^ — > /i^ + ; 

n + n + 7r+^ /i^ + ; 
—i-p + n + Tr'^^j + j 



(4) 



which occur in approximately equal ratios and with near 
unit total probability. The corresponding p -I- p {n-\-n) — > 
TT* , 7r° processes do not involve a relative drift velocity 
(as do the p -I- n), and are thus less probable. Processes 
leading to multiple baryons are also suppressed due to the 
higher threshold, and for simplicity we restrict ourselves 
to the above p + n processes. 



III. 10 GEV NEUTRINOS AND 7-RAYS 



The total number of neutrons carried by the fireball is 



r]7npC^ 
>53 



(5) 



The comoving optical depth r' ~ ripCrr/F oc a/{rT) has 
the same dependence for pion formation and photon scat- 
tering, but cr^ <C (Tt (Thomson cross section), so the pi- 
onosphere r-j^ occurs below the 7-photosphere r^. The 7- 
rays in equation (^) can only escape from a skin depth be- 
low the 7-sphere in the essentially laminar flow with prob- 
ability <, T^ir.^) - r^/r.^ - ((7^/aT)(l+?/7)^ - 1/25, 
for 77 ^ j/tt- Each n leads to '--^ 1 photon of CM en- 
ergy e'^ ~ 70 MeV and observer energy centered broadly 



around e-y ~ 70F„//(1 + z)MeV - 10 GeV. Using a 
proper distance Dp — 2.8 x lO'^^hgr^ll — + z] cm 

with a Hubble constant /igs = i?o/65Km/s/Mpc, the 
number fluence at Earth is ~ N„P^/ATrDp ~ 10~^ 
(2) cm^^. This is below the sensitivity of the ^ 200 cm^ 



area EGRET detector on the Compton Gamma Ray Ob- 
servatory (e.g.l3), but ipr rare nearby bursts it may be 
detectable by GLASTE 
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The neutrinos originate at the pionospheric radius 
<C where t^^ 1. In this region the stable charged 
products and 7-rays from the reactions (|^) remain in the 
fireball, and each n leads on average to one v and one 
P. We list below the average neutrino energies for pions 
and muons decaying at rest. The neutrinos from muon 
decay have a continuum spectrum. Also, the energies are 
Doppler broadened by Vrei/c 0.5. 

e!-, ~ 30MeV , e' ~ SOMeV from 7r± 

e'^^ ~ SOMeV , e'p^ ~ SOMeV from /i+ , (6) 

ep_^ ~ 30MeV , e'^^ ~ 50MeV from ^- 

The relevant cross section for detection averaged over v 
and V ^^<J^ ^ 0.5 x 10~'^*(e/GeV) cm^ at the observed en- 
ergy en The observer frame energy is e = e'Q;r„^/(l-|-z), 
where a ~ 1 near threshold. For the CM vv production 
energies of equation (|^), the average v + v CM energy 
per neutron is e' ~ 100 MeV. Taking a ~ 1 the observer 
V -\- V energy per neutron is e ~ 0.ir„//(l + z) GeV, 
and the effective detection cross section per neutron is 
CT^p 5 X 10~''°r„/(l + z)"^ cm^. Multiplying by a burst 
rate within a Hubble distance oiTZb ^ lO^T^-bs/year, for a 
1 kni'^ detector containing Nt ^ \Q^^Ntz9 target protons, 
the rate R^p = (A^t/47rDp)7?.6iV„CT^p is 




events in the detector in coincidence with GRB electro- 
magnetic flashes. The energies of the events are 

e,^p^ - lOGeV , e,,^, ~ 5GeV , (8) 

which scale (x EliX\'\~J'' (2/[l + S]Y'\2/[l + z]) 

Subsequent to decoupling and n,p collisions, each neu- 
tron decay n ^ p + e~ -\- leads to an additional of 
CM energy e^^ ^ ^ 0.8 MeV, which boosted in the ob- 
server frame by r„//(l -I- z) is ^ 12QMeV. The cross 
section is Cp^ ~ 2 x 10"^" cm^ and the expected rate in 
a km'^ detector is less than one event per year. 

IV. MEV 7-RAYS 

The non-thermal MeV "f-iB^s are thought to be pro- 
duced in collisionless shocks,liil which occur at a radius 
rsh > rj > Ttt, after the bulk Lorentz factor has satu- 
rated to its asymptotic value. For an n,p outflow, shocks 
can occur in the original p, as well as in the n component 
after the latter have decayed, and this can influence the 
external shock light curves.ll2l A separate and important 
consequence of neutron decay is that it should also affect 



the internal shock gamma-ray light-curves. In the proton 
component internal shocks occur at Tsh ~ ct^Fp^, where 
ty is the variability timescale, and Tpf is the asymptotic 
proton Lorentz factor. From energy conservation, for 
77 > ,7, this is rp/~77(l+0[l-(e/[l+?])(6/7)(r7,/r7)4/3], 
and taking into account photon drag one can show that 
an upper limit is Tpf,„iax ^ 8.3 x lO'^ElJ^'^Q^^'^r^^^'^ 
(l + 8^/7)-i/4. 

The 7-rays start to arrive at an observer time t ^ ^ 
10~'^t.u_3 s, lasting for a time t^, (where 10"'^ ^ ^ lO'^ 
s). For the n component, r^h ctmi-nXl^f, with r„/ 
from equation (||) and i,„i„ ~ min[ty,t'^V~^], where t'^ ~ 
lO'^ s is the comoving frame neutron decay time. Taking 
^ ~ 1 in the estimates below, for 20 ^ 77 ^ 77^ ~ 400 
the neutrons decay and shock beyond the proton shock 
for any < 10'^r7~^, at observer times i„ = [50s, 3s], 
while for f] ?2 1],^ ^ 400 the neutrons decay and shock 
beyond the proton shock for any <5 3(r;^/r;)^/'^ s. The 
typical observed duration of the decay, including the blue 
shift due to the bulk motion towards the observer, is 
tn ^ 10^/F„/, where F„/ ~ 77 for 77 < 77^ and F„/ = 
(3/4)77^(77/7;^)^^/'^ for 77 > 7;^. Thus t„ decreases from 
approximately 50 s to 3 s for 20 ^ 77 ^ 77^ ~ 400, and 
then slowly increases again as t„ ^ 3{ri/r]jry/^ for rj ^ 
tItt ~ 400, with both 777r and t„ scaling cx [(1 + S,)/2]^^/'^. 

The number of neutron decays is oc 1 — exp(— t/t„), 
so the envelope of the neutron-related light curve is the 
mirror image of a "fred" (fast rise - exponential decay) , 
i.e. an "anti-fred" (or generally, slow rise - fast decay). 
In general, photon emission starts at from the proton- 
related component, which lasts a time with an arbi- 
trary shape envelope, modulated by spikes of minimum 
duration ty < t^, depending on the chaotic behavior of 
the central engine producing the outflow. The neutron- 
related component starts at a later time i„ > ty, and 
has an anti-fred shaped envelope modulated by spikes 
of ty and a total duration tyj. If > t„, the anti-fred 
component would be hard to distinguish because of the 
superposition of the ongoing p and 77, components. How- 
ever, for short bursts with < t„ ~ 3 s, the p and n 
components are separated: first there is a pulse of du- 
ration ty, with a random envelope, followed after a time 
~ tn by a pulse with an anti-fred envelope of duration 
tn, and characteristic photon energy softer than the pre- 
vious by e„/ep ~ ty,/tn (which if small could be below 
the BATSE band, but may be detectable with the Swift 
satelliteE3) . The latter pulse is a signature for neutron 
decay in the burst. 

V. DISCUSSION 

For characteristic parameters, GRB outflows produce 
5-10 GeV v^j.D^ and i/gi/e from internal inelastic p, ?t. col- 
lisions that create pions. The vv energy output ^ 
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5 X 105ii;53(2e/[l + C])(2/[1 + z])(?7^/?7)4/3 ergs depends 
on the total energy E of the GRB and on the neu- 
tron fraction ^ as weh as on the dimensionless entropy 
r]. For a km^ detector, approximately 5 to 10 neutrino 
events above 10 GeV are predicted per year, for a neu- 
tron/proton ratio ^ = 1. These events will be coinci- 
dent with GRB electromagnetic flashes in direction and 
in time (to an accuracy of ~ 10 s), which can enable their 
separation from the atmospheric neutrino background. 
Underground-water detectors of the type being planned 
by BAIKAL jia NESTOBid ANTARESta and the Antar- 
tic detector ICECUBEiij could potentially detect these 
relatively low energy neutrino events if a sufficiently high 
density of phototubes were used. About 80% of these 
neutrinos are and (in approximately equal num- 
bers) and the remainder are i^e and Uf,. These 5-10 Gev 
lyp are followed by ^ 120 MeV De from neutron decay, 
but the event rate from neutron-decay neutrinos is very 
low. The higher energy neutrinos are produced for neu- 
tron/proton ratios ^ > when the dimensionless entropy 
77 = L/Mc2 exceeds t]^ ~ 4 x lO^Lg/V^^^^"* (2/[l + 
and are accompanied by ~ 10 GeV photons which may 
be detectable in low redshift cases with GLAST.eI For a 
typical GRB at redshift z ^ 1 the number fluences in 10 
GeV neutrinos are N{De + Ve) ^ 0.5A^(P^ + v^) ^ 10~^ 
cm~^, and one order of magnitude less for GeV photons. 

In all bursts where ^ > the lower energy (~ 120 MeV) 
neutrinos are produced, and neutron decay occurs on an 

observer timescale i„ ~ ?,L~^''^rllj'^[{l+C)l'A'^''^{r]lri^Y'^ 
s. For outflows of duration tw, these decays will be as- 
sociated with MeV electromagnetic pulses of duration 
min[i„,tu,], which are additional to the MeV pulses ex- 
pected from shocks in the original proton component. 
For short bursts with <5 3 s, the proton electromag- 
netic pulse appears first and is separated from a subse- 
quent neutron electromagnetic pulse, the latter having 
a slow rise-fast decay envelope and a softer spectmm, 
which may be detectable with the Swift satelliteJld A 
systematic study of the time histories of GRB emission 
would be useful to search for evidence of delayed pulses 
that might be caused by neutron decay. 

The detection of 5-10 GeV vv in coincidence with GRB 
photon flashes will not be easy, but would provide unique 
astrophysical information. Gonstraints on the neutron 
fraction could provide information about the progenitor 
stellar system giving rise to GRB. For instance, core col- 
lapse of massive stars would lead to an outflow from an 
_Fe-rich core with ^ ~ 2/3 — 1, while neutron star merg- 
ers would imply ^ > 1. Photodissociation during collapse 
or merger, as well as n,p decoupling and inelastic colli- 
sions, would both drive ^ toward unity, although this 
equalization process is likely to remain incomplete. For 
low 7? ^ ?77r, inelastic collisions are not expected and the 
5-10 GeV vv are absent, producing only the harder to 
detect ^ 100 MeV Ve from neutron decay. An initially 



non-baryonic outflow of, e.g. e and magnetic fields, 
would acquire a baryonic load by entrainment from the 
progenitor environment, with ^ <C 1 from massive stellar 
envelopes, but ^ ^ 1 for, e.g., compact mergers. Thus, 
lower values of ^, leading to lower ratios of 5-10 GeV vv 
and a lower ratio of neutron decay MeV photons to to- 
tal fluences would be expected from massive progenitors 
than from compact mergers. 
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